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ABSTRACT 

Anisotropic gravitational radiation from a coalescing black hole binary can impart a recoil velocity 
of up to several hundred km/s to the remnant black hole. We examine the effects of recoiling massive 
black holes on their host stellar bulges, both for holes that escape their host and those that return to 
the galactic center via dynamical friction. We show that removal of a black hole via radiation recoil 
generally results in a rapidly-formed central core in the stellar system, with the effect being largest 
when the hole stays bound to the bulge and the recoil velocity is comparable to the bulge velocity 
dispersion. Black hole recoil therefore provides a mechanism for producing cores in some early type 
galaxies, but it is expected to be most efficient in faint ellipticals that are known to have steep density 
profiles. We argue that these results may hint at a significant role for gas in facilitating the coalescence 
of binary black holes in faint (power-law) early-type galaxies. 
Subject headings: 



1. INTRODUCTION 

When galaxies with central black holes merge, dynam- 
ical friction will drive the holes toward the center of 
the remnant, creating a black hole binary. If the bi- 
nary separation becomes small enough that gravita- 
tional wave emission is significant, the binary will rapidly 
coalesce. The primary mechanism for reducing a;, is 
three-body interactions: stars with orbits passing close 
to the binary are scattered to high velocities, allowing the 
binar y's orbit to decay ijBegelman. Blandford. &: ReesI 
Il980|) . The timescale for removing stars from low an- 
gular momentum orbits (the "loss cone") via gravita- 
tional slingshot is much shorter than the (coUisional) 
timesca le for repopulating the loss cone in a spherical 
galaxy iMakino fc Funatol200l. The lat ter is many Gyr 
or longer l|Milosavlievic fc Merrittll20d3^ ■ implying that 
the decay of the binary would stall before gravitational 
waves become important. 

Recent work has found that a number of pro- 
cesses can increase the hardening rate: Brown- 
ian motion of the bin arv l| Quinlan fc HernauistI Il997t 
[Chattcriec, Hernauist, fc Loebl l2003j) . triaxialitv of the 
stellar bulge ( Yu 2002) , interaction with a massive accre- 
tion disk (Armitage & Natara.j an 2002), and re-ejection 
of p reviously ejected stars l)Milosavlievic fc Merritti 
l2ml all may help the binary coalesce. Coalescing black 
hole binaries radiate gravitational waves that carry away 
energy and both angular and linear momentum. The 
remnant black hole then receive s a recoil velocity fkirk in 
the range of 100 to 500 km/s ijFavata. Hughes, fc Hold 
1200^ . The black hole can therefore either leave the 
host galaxy (if iikick > Wosc) or oscillate around the cen- 
ter of the galaxy, its orbit decaying due to dynamical 
friction. Radiation recoil has a number of astrophysi- 
cal implications, from constraining the growth of high- 
redshift quasars to the possibility of an intergalactic pop- 
ulation of massive black holes (Madau fc Ouataert 2004; 
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iMerritt et 311120041: lHaimanll2004|) . In this paper we use 
numerical experiments to study the reaction of stellar 
bulges to co alescing black holes at their centers (also see 
IMerritt et a,1...20041 



2. SIMULATIONS 

We perform N-body simulations to follow the evo- 
lution of a stellar bulge containing a central super- 
massive black hole that has been given a recoil veloc- 
ity immediately following a binary coalescence. The 
purely gravity) simulations are performed using GADGET 
Sprin gcl. Yoshida. & White 2001), a publicly available 
N-body tree code. We generate initial conditions us- 
ing the equilibrium distribution functio n for a spherical 
stella r system with a central black hole l)Tremaine et alJ 
Il994j) . Since dark matter contributes a small fraction of 
the mass in the inner tens of parsecs of a typical ellipti- 
cal galaxy, we ignore its contribution in our simulations. 
We use the iHernguist (1990) model for the stellar sys- 
tem; this model is attractive because of its simple form 
and the fact that in projection it closely resembles the de 
Vaucouleurs R^/^ surface brightness law. Furthermore, 
energy transfer from orbital decay of a black hole binary 
to surrounding stars is expected to transform dense cusps 
associated with the gr owth of a black hole to profile s no 
steeper than p oc (jMilosavlievic fc Merrittll2001j) . so 
an cusp serves as an upper limit for the expected stel- 
lar distribution around a coalescing black hole binary. 
The density profile of the Hernquist model is given by 



p{r) 



M» a 



2T:a^ r (1 + r/a)^ 



with total mass A/* and scale radius a, related to the half- 
mass radius by ri/2 = (1 -I- V2)a. The escape velocity 
from the center of the system and the dynamical time and 
circular velocity at the scale radius are given by: Uosc — 
(2GM,/a)i/2, tdyn = tdynir = a) = {Sn^ayGM^y/^, 
and Vch-c = Wcirc(?' — a) = (GM^/AaY^^ . We consider 
only spherically symmetric, isotropic models, in which 
case the phase-space distribution function depends on 
energy alone and can be calculated from Eddington's for- 
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mula l)Binnev fc Tremain3ll987j) : 

Here ip is the negative of the gravitational potential, 
, GM, , GMbh 
r + a r 

and E is the binding energy per unit mass {E > 0). We 
note that without a central black hole {Mbh = 0), the 
only accessible energies are < -E < GM^,/a and f{E) 
is analytic. Adding a black hole allows particles to have 
E > GM^ja and j(E) must be computed numerically. 
We initialize particle positions from the mass profile, and 
then obtain each particle's energy from the calculated 
distribution function: P{E\r) oc f{E)y/^{r) — E. 

In the presence of a central point mass, f(E) is 
not a monotonic function of energy fsee iTremaine et al.l 
119941 for plots of f{E) for several values of Mbh/M^). 
Stability to radial perturbations must therefore be de- 
termined numerically. We have tested that our sim- 
ulations (without black hole recoil) indeed maintain 
equilibrium for many dynamical times. Such sys- 
tems are expected to evolve on the t wo-body re- 
laxation timescale (jSpitzer &: H art' '1971'), forming a 
p (X r^'^l'^ cusp on small scale s (Bahcall & Wolf 1976; 
iPreto. Merrit. fc Spurzem 20041 , but the relaxation time 
is of order 10^° years or longer in a typical elliptical 
galaxy, much longer than the timescales relevant to this 
paper. 

We set the black hole to stellar bulge mass ratio to be 
MbhI^'^* — 1/300 and simulate cases with z^kick both 
above and below Ucsc for comparison. We use iV = 10^ 
equal-mass particles to represent the stellar bulge, and a 
force softening of e = 0.0176a that allows us to resolve 
down to ~ 0.2rft,, where r-^ is the sphere of influence of 
the black hole given by M(< r^) = 2Mbh (a definition 
equivalent to the standard r^j = GMbh /o'^ for a singu- 
lar isothermal sphere) . is typically a small fraction of 
the scale radius of the stellar system, e.g., = 0.089 a 
for Mbh/M^ = 1/300. Since gravity is the only physics 
in the simulations, our results can be interpreted at dif- 
ferent length and mass scales by rescaling: for a — > Aa 
and M /3M, take v (P/X)^^^ u and i -> (A^//?)^/^ t. 
For reference, a model with stellar mass 4 x 1O*M0 and a 
scale radius of a = 42.5 pc has tdyn=1.15 Myr, Wcirc = 100 
km/s, Vesc — 283 km/s, and a projected central veloc- 
ity dispersion at one-eighth the effective radius Re oi 
ap{RE/S) « 70 km/s. 

3. STELLAR DENSITY PROFILES 

A consequence of gravitational radiation recoil is evo- 
lution of the density profile p{r) of the stellar system. 
Fig. n] compares p{r) at three epochs soon after the re- 
coil of the black hole. Within O.l^dyn after the black 
hole leaves the central region of the bulge, p{r) flattens 
substantially. A core of ~ rh forms and remains for the 
entire length of the simulation. 

Three primary mechanisms could contribute to the 
rapid flattening of the inner stellar density profile: (1) 
departure of stars bound to the black hole; (2) energy 
deposition due to dynamical friction on the black hole; 
and (3) re-equilibration of the stellar system due to the 
departure of the black hole. We discuss each one in turn. 
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Fig. 1. — Early evolution of the stellar density profile (in units 
of M./a^) after the black hole (with Mbh/M* = 1/300) re- 
coils (at t = 0). Three simulations are compared: recoil veloci- 
ties above (dashed) and below (dotted) the bulge escape velocity 
VcBc = 2.83ticirci S'lid test model in which the hole is removed 
instantly at t = (diamond). All three runs produce a density 
core inside ~ within O.OStfdyn after the recoil; p evolves little 
afterwards, so for clarity we offset p at 0.43 and 0.87 fdyn by a 
factor of 3 and 8. For reference, if M* = 4 X 10* Mq and a = 42.5 
pc, tdyn=l-15 Myr and r,^ = 0.089 a = 3.78 pc. 



In process 1 , some stars move with the black hole as it 
recoils from the galactic nucleus. The mass of stars Mt, 
bound to the black hole is a strongly decreasing function 
of the recoil velocity: for Ukick/fcirc = (0,0.75,1,2,3), 
our simulations yield Mb/ Mbh = (0.769, 0.394, 0.249, 
0.0157, 0.0012). The removal of stars bound to the re- 
coiling black hole is therefore too small to account for 
the mass deficit of ^ 2Mbh seen in the simulations. 

Dynamical heating due to energy transferred to stars 
from the decay of the black hole's orbit (process 2) is 
relatively slow. The heating timescale is ~ \E/E\, which 
is about 5.7tdyn (in the Wkick — 2fcirc run) given that 
the change in the black hole energy from start to apoc- 
enter is AE/Et=Q ~ —0.076 and that the time required 
for the black hole to reach apocenter is « 0.43 tdyn- In 
comparison. Fig. ^ shows core development over a much 
shorter timescale (< 0.1 tdyn), indicating that dynamical 
friction heating is not the dominant process in the initial 
core formation. Its longer term effect will be discussed 
below. 

To study the effects due to process 3, we first note that 
the timescale for the stellar system to re-adjust after the 
black hole ejection is the crossing time of the core, which 
is « 0.075 tdyn for a core of « 0.14a shown in Fig. ^ 
This timescale is comparable to the core formation time 
of < 0.087 tdyn; the stellar system therefore has suffi- 
cient time to respond to the black hole kick by producing 
the observed fiattening of p(r). We test the hypothesis 
that process 3 dominates the early evolution of the stellar 
system by performing a simulation in which we remove 
the black hole instantaneously from an initial equilib- 
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Fig. 2. — Evolution of the stellar density profile at 0.87, 1.74, 
and 3.48 t^yn foi" the same runs as in Fig. 1. For i^kick = 3i>circ > 
Dcsc, the density profile remains constant once the hole leaves the 
galactic nucleus, while for Dkick = 2 ticirc < ^csc there is both early 
evolution as the stellar bulge dynamically adjusts to the absence 
of the hole and later evolution due to dynamical friction heating 
as the black hole returns to the center of the stellar system. 
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Fig. 3. — Comparison of the ii^ick = 2vcirc < '^csc simula- 
tion (dotted curves) with a calculation in which the hole is in- 
stantaneously removed at t = and re-introduced at 0.87 t^y^^. 
The results highlight the role of dynamical friction: the static re- 
introduction of a black hole to the core causes the re-building of 
a central cusp by t = 1.74ttjyn, while the true simulation shows 
core-formation by dynamical friction heating as the hole returns to 
the galactic center. 



rium configuration and then allow the system to evolve 
to its new equilibrium (also see Merritt fit al. 2004). This 
test eliminates evolution due to gravitational interactions 
between the black hole and the stars and mimics the ef- 
fect of Wkick ^ Vcsc- Fig. n compares the stellar p{r) 
for this model with the Ukick = 2 Wcirc and 3 Ucirc simu- 
lations. This figure shows that at early times, all three 
simulations behave quite similarly; in fact, p(r) is virtu- 
ally indistinguishable at both 0.087 and 0.43 idyn- This 
supports the hypothesis that dynamical adjustment to 
the absence of the black hole is responsible for the initial 
core formation. 

The longer-term evolution of the stellar density dis- 
tribution is shown in Fig. [21 For the run where the 
black hole escapes (wkick = 3wcirc), pi^) is essentially un- 
changed from 0.43 idyn up to the end of our simulation at 
3.48 idyn, at which point the black hole is more than 13a 
from the center of the stellar system. The evolution is 
qualitatively different, however, if the black hole returns. 
For i^kick = 2ucirc, the core density at 0.87 tdyn is higher 
than that with z^kick > "^esc- This effect is due to the 
presence of the black hole (and the stars bound to it) on 
its return path through the center of the system. Subse- 
quently, however, the core becomes less dense relative to 
the Ukick = 3ucirc run, at both 1.74 and 3.48 tdyn- From 
3.48 tdyn onwards, p{r) of the i^kick = 2wcirc run is essen- 
tially unchanged and the black hole is nearly stationary 
at the center of the stellar potential. In both simulations, 
the velocity dispersion tensor remains isotropic and the 
bulge remains mostly spherical. 

To explore the origin of the additional flattening in 
the Ukick — 2 Wcirc Calculation, we perform a control sim- 
ulation in which we re-introduce the black hole into the 
"instantaneous removal" simulation at 0.87 tdyn- In this 



case, dynamical friction is irrelevant (as in the instan- 
taneous removal simulation) but the stellar system re- 
adjusts because of the black hole. Fig. |21 shows that 
the equilibrium profile attained after re-addition of the 
black hole is significantly steeper than that found in the 
Vkick = 2wcirc simulation at t = 1.74 idyn- This compari- 
son test shows that the additional flattening seen in the 
i^kick — 2 Wcirc run is due to the heating of the stellar sys- 
tem by the black hole as dynamical friction returns it to 
the center of the galaxy. Even though this input of en- 
ergy is less than 1/1000 of the total energy of the initial 
stellar system (for Wkick — 2wcirc)j it is enough to ensure 
that the new equilibrium differs substantially from the 
original at small radii. 

The results of this section show that for recoil velocities 
I'kick I'circ, a stccp density cusp is difficult to maintain 
in a stellar system containing a coalescing black hole bi- 
nary (simulations with ■ykick ^ 'i'circ show little change in 
the density profile of the bulge). We now turn to the 
implications of these results. 

4. DISCUSSION AND IMPLICATIONS 

In a hierarchical cosmology, mergers of galaxies with 
central black holes will lead to the formation of black hole 
binaries. This process occurs on timescales significantly 
shorter than the Hubble time only in the case of ma- 
jor mergers, where the ma ss ratio of the two galaxies is 
withi n a factor of about 3 l)Volonteri. Haardt. fc MadaiJ 
|2003|) . If the timescale for coalescence of the binary black 
holes is shorter than the typical time between galaxy 
mergers, it is likely that gravitational radiation recoil 
will have interesting observational consequences, most 
notably for early type galaxies. 
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Observations of early-type galaxies show that 
their density profiles fall into tw o dist inct classes 
ijFerrarese et alJ Il994t iLauer et Power 
law galaxies, which tend to be faint, have surface 
brightness profiles that are steep and do not seem 
to exhibit a significant break at small radii. Core 
galaxies are generally brighter and have a discernible 
break in their surface brightness profiles. Faber et 
al. (1997) discuss several possibilities for the origin 
of core galaxies. In particular, they and others have 
argued that the formation of a black hole binary can 
produce a core in a stellar system, both by gravitational 
slingshot of stars during the decay of the bina ry and 
by Brownian motion of the binary (| Ebisuzaki_et_aLj 
1991; Makino & Ebisuzaki 1996; Ouinlan & Hcrnauist 
1997: Faber et a l. 1997; Milosavlicvic_&_Merritt 
20011: iRavind ranath. Ho & FihoDcni^l 120021 

Milosavlievic et alJ i2002,) . This process of central 
density reduction due to loss cone depletion should be 
ubiquitous in galaxies with binary black holes, reducing 
their central cusps to p cx or shallower. 

We have demonstrated in this paper that in addition 
to these mechanisms for core formation, gravitational re- 
coil in the predicted velocity range of 100 to 500 km/s 
generically results in the formation a core in the density 
profile of a stellar system, which produces a core in the 
surface brightness analogous to that seen in bright ellipti- 
cals. (We note that the late-time volume density profiles 
in Figs. lllHl also yield surface densities that are nearly flat 
inside « 1 — 2 r/,, but do not affect the projected aper- 
ture velocity dispersions noticeably.) For recoil velocities 
^kick > ^^esc = 2.83uciic (for thc Hcrnquist model), which 
are rare for Milky Way size bulges or larger but plausi- 
ble for smaller galaxies, we expect the black hole recoil to 
produce stellar cores that are largely independent of Wkick 
since stellar re-equilibration is the dominant process. For 

t'kick ^ "^^esc — 2.83l'circ; 

however, recoils of a given Wkick 
have a larger effect on the inner stellar profile in smaller 
galaxies because Wkick/wcirc is larger. In addition, black 
holes that do not escape their host galaxy tend to create 
larger cores via dynamical friction than black holes that 
do escape. In the context of these results, it is thus the 
existence of power law galaxies that must be explained. 
Indeed, one might naively have predicted that smaller. 



fainter galaxies would have proportionally larger cores 
(in units of r^) because of gravitational recoil, contrary 
to what is observed. 

One resolution of this difficulty is that black holes are 
rare in low-mass galaxies. The dynamical constraints 
on black holes in such galaxies are indeed rather poor, 
though there are suggestive hints of substantial popu- 
lation of ~ 1 0^ - lO^MfT) bl ack holes in dwarf Scyfcrt 
galaxies (e.g. iFilitrnenko fc H o 2003; Earth ct al. 2004 
iGreene fc Holl2004|) . There are also plausible mecha- 
nisms that could rebuild a cusp. A natural explana- 
tion for reforming a power-law surface brightness profile 
is star formation accompanying a dissipative merger of 
gas-rich galaxies. This scenario requires the black hole 
binary to coalesce - and the recoiling hole to return to 
thc galactic nucleus - before the starburst is completed. 
If stellar dynamical processes alone are responsible for 
hardening the binary, the timescale is typically > 10^ 
years even in the case of highly triaxial galaxies l)YiJ 
|20Q2)- This is comparable to, or somewhat longer than, 
the e xpected duration of a merger-triggered starburst 
('e.g.. lMihos"fc Hcrnauist 199^. An intriguing possibil- 
ity is that the presence of substantial amounts of gas can 
greatly reduce the time for binary black holes to coalesce. 
If this is the case, as argued by, e.g., Armitage & Natara- 
jan (2002) and Escala ct al. (2004), then it is plausible 
that a black hole binary can coalesce before the end of a 
starburst and a stellar cusp can be re-formed, thus help- 
ing explain the existence of power-law early-type galax- 
ies that are built hierarchically. We note that, given this 
scenario, in future work it would be interesting to con- 
sider the dynamics of recoiling holes in the presence of a 
surrounding accretion disk. 
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